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^ • Abstract 

cn ' We investigate the stability against inhomogeneous perturbations and the appear- 
ance of ghost modes in Gauss-Bonnet gravitational theories with a non-minimally 

psj I coupled scalar field, which can be regarded as either the dilaton or a compactifica- 

K^ ' tion modulus in the context of string theory. Through cosmological linear pertur- 

^—. • bations we extract four no-ghost and two sub-luminal constraint equations, written 

fSJ ! in terms of background quantities, which must be satisfied for consistency. We also 

^T I argue that, for a general action with quadratic Riemann invariants, homogeneous 

vQ ' and inhomogeneous perturbations are, in general, inequivalent, and that attractors 

f^ . in the phase space can have ghosts. These results are then generalized to a two-field 

^ I configuration. Single-field models as candidates for dark energy are explored numer- 

^ • ically and severe bounds on the parameter space of initial conditions are placed. A 

^H. number of cases proposed in the literature are tested and most of them are found 

r~| ! to be unstable or observationally unviable. 
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1 Introduction 



General Relativity (GR) is a very accurate and successful theory of classical 
gravity which would be desirable to embed in a more fundamental theory. 
The formulation of a quantum theory of gravity has proven to be a very 
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difficult task and, so far, string theory is the only framework (together with 
loop quantum gravity) within which gravity can be incorporated at quantum 
level in a consistent way. Therefore it is worth exploring low-energy aspects of 
gravity considering, at the same time, leading-order string corrections. This 
leads us to the study of low-energy actions that contain a Gauss-Bonnet (GB) 
term, which is a particular combination of quadratic Riemann invariants. 

One of the most striking peculiarities of string theory is the fact that all 
couplings, including those to GB terms, are given dynamically, i.e. as vacuum 
expectation values of fields (called, in general, moduli). This is a feature which 
makes the study of the evolution of the system of gravity plus moduli fields 
particularly interesting, as done for a Friedmann-Robertson-Walker (FRW) 
background in [1]. In four dimensions, the GB term is topological and does not 
contribute to the dynamics. Hence it is necessary for it to have a nonconstant, 
moduli-dependent coupling if we want it to play a nontrivial role, unless one 
appeals to a Lorentz-violating configuration such as the braneworld. 

In this article we deal with the cosmological study of GB gravity and focus on 
a usually overlooked aspect, namely the possibility of having ghosts and other 
quantum instabilities in the model. A ghost is, by definition, a field whose 
kinetic term in the action is unbounded from below (roughly speaking, it has 
the 'wrong sign'), which implies both a macroscopic instability which, if not 
healed, would lead to a breakdown of the theory, and a violation of unitarity. 

It is of utmost importance to realise that the presence and nature of ghost 
modes for a covariant Lagrangian depends on which classical background the 
theory is assumed to live in. For instance, in string theory the action is natu- 
rally expanded around a Minkowski target space. Ghosts in minimally coupled 
GB and higher-derivative gravities in constant curvature backgrounds were 
studied in [2,3,4,5,6]. Those approaches include both the Minkowski (ffist ex- 
plored in [7]) and de Sitter (dS) cases, and infer the absence of spin-2 ghosts in 
actions whose higher-order term is a generic function of the GB combination. 

One might be tempted to assume that, since de Sitter is a ghost-free vacuum 
of the GB theory, any cosmological model with such Lagrangian does not suf- 
fer from quantum instabilities. However, the FRW curvature invariant R is 
not constant and the methods used in the above-mentioned papers cannot be 
applied any longer. Some works have pioneered this issue by studying cosmo- 
logical (i.e. FRW) perturbations for non-minimally coupled actions [8,9,10], 
albeit the ghost problem is not explicitly addressed there. 

In analogy to the situation when one considers scalar-tensor theories of gravity, 
where the absence of spin-2 ghosts require constraints on the vacuum expecta- 
tion value of the scalar (see [11,12]), we have extracted a number of constraints 
on the field-dependent couplings of a GB model with single scalar field in or- 



der to avoid the appearance of instabilities. These are obtained by computing 
the gravitational perturbations about a FRW background and studying the 
scalar, vector and tensor contributions separately, as done also in [8] and sub- 
sequent papers. We also imposed that perturbation modes do not propagate 
faster than light. 

We have then applied the obtained results to different single-field models in 
literature [1,13,14], whose stability was studied only in phase space, that is, 
at the classical level (see also [15]). The no-ghost and suh-luminal constraints 
are here regarded as strong selection rules in the space of parameters of the 
theory, and we provide an algorithm necessary (but not sufficient) to guarantee 
quantum stability and obtain a suitable acceleration today. 

We have found that, for general initial conditions and during cosmologically 
relevant periods, ghost modes arise in most of the string-inspired models stud- 
ied so far. ^ By fine tuning the GB coupling constant /?, however, it is possible 
to get viable models free of instabilities and reproducing the observed accel- 
eration today. In particular, j3 must be so small that the GB term is actually 
unobservable at late times. 

If one regards the GB action as an effective one, there might be higher-order 
curvature corrections which are subdominant relative to the leading GB term 
at late times. Hence, if the latter gives rise to an unstable theory, the higher- 
order ones will do so as well (the no-ghost conditions would still be dominated 
by linear and quadratic terms, with no appreciable impact from higher powers 
in /9 <^ 1). Conversely, at high energies or very early time, the effective theory 
might have received higher-order or nonperturbative modifications which are 
presently out of control in the sense that, from the study of only the GB term, 
we cannot know whether they are ghost free or how the universe evolved due 
to their action. Then one should take this other Lagrangian, whatever it is, 
and study it from the very beginning as a separate case. The key point is that 
the GB action below is possibly inadequate, by itself, for describing cosmology. 

The presence of attractor solutions in phase space does not guarantee by itself 
quantum stability, and vice versa. In the context of an action with a general 
quadratic combination of Riemann invariants, we provide general arguments 
and several concrete examples showing that stability against homogeneous 
(i.e. classical in phase-space sense) perturbations is unrelated to that against 
inhomogeneous (i.e. cosmological, quantum) ones. ^ 



^ A parallel study of modified gravity models can be found in [16]. These insta- 
bilities do not spare even more complicated cosmological models with an arbitrary 
function of the GB combination [17]. 

^ Throughout the paper we refer to cosmological perturbations as 'quantum' but 
in fact they are still (semi)classical, since the perturbed equations do not incorpo- 
rate loop interactions between particles. However, inhomogeneous perturbations are 



If probability is not conserved in scattering processes, as when ghosts are 
involved, wild particle creation can occur in a time interval of cosmological 
length. For this reason, we define as physically viable those solutions that do 
not have ghost modes at any time during the evolution of the universe. This 
attitude differs somewhat from other interpretations in literature [9] , which we 
shall briefly discuss in Section 5.4. Although the presence of unstable modes 
may suggest that the adopted background (FRW) is pathological, whereas the 
theory per se (as covariantly formulated) is not, nonetheless this background 
is the best known model reproducing natural, observed phenomena on large 
scales, theoretically consistent with the inflationary scenario. Hence, we will 
consider solutions to be physically unviable if they lead to ghost modes and/or 
causality violation on FRW. Therefore, here we do not address the perhaps 
more abstract problem of the 'vacuum structure' of these theories. In our 
opinion, inconsistencies in FRW are sufficient to discriminate among models. 

In [18,19] an estimate for the decay rate of minimally coupled ghost particles 
into photons via gravitational interaction was calculated. While theories with 
a Lorentz invariant cut-off are experimentally forbidden, there remains the 
possibility of a Lorentz-violating cut-off Auv ^ 3 MeV. In our framework, 
ghost modes are non-minimally coupled to gravity, but this can only tighten 
the bound of [19] (see that reference for details). Moreover, the characteristic 
energy scale of the graviton-modulus interaction vertex in modulus-driven 
cosmologies is that of string theory, hence much higher than Auy. Then each 
case in which we encounter a ghost mode is automatically excluded. Also, we 
shall use a very compressed time variable r which highlights only features 
extended through cosmologically long time intervals. The simple fact that one 
can see the ghosts in r-plots means that there is no hope for the effect of such 
instabilities to be unobservable. 

The paper is organized as follows. In Section 2 we write down the covariant 
and cosmological equations of motion describing theories with a Gauss-Bonnet 
term and one non-minimally coupled scalar field. In Section 3 we present the 
general structure of the solutions to the cosmological equations outlined in 
the previous section, describing some numerical issues associated to these so- 
lutions. In Section 4 we study the tensor, vector and scalar perturbations 
about a FRW background and derive four conditions for the absence of quan- 
tum instabilities in Gauss-Bonnet cosmologies, plus two more conditions on 
the speed of propagation of tensor and scalar modes. These conditions are ap- 
plied to some string motivated models existing in the literature in Section 5. 
The problem of suitable initial conditions, giving stable and experimentally 
viable solutions, is considered in Section 6; the theory is here regarded as a 



closer to the quantum picture than homogeneous ones, in the sense that they iden- 
tify particle modes and their tree-level dynamics. As said above, ghost instabilities 
in general extend beyond the classical level anyway. 



model which describes the present acceleration of the universe as measured by 
supernovae and large-scale structure observations. The influence of other ad- 
ditional scalar degrees of freedom coupled to the GB combination is discussed 
in Section 7. Section 8 is devoted to the discussion of our main results. 



2 Setup 

We are going to study gravity models with a Gauss-Bonnet term, that we 
write in the standard form 

Cgb = Q-^P + R\ (1) 

where the Riemann invariants 

P = RfiuR^^ 1 Q = RjiuXkR^^ '^ , (2) 

are built out of the Riemann tensor, Ricci tensor and Ricci scalar which are 
defined as 

pA -- fj -pA o -pA _|_ -per pA -pa pA /o\ 
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Rp-v = R uxv, R = R^ivg'"' , (4) 



where T^^ are the Christoffel symbols. Greek indices run from to 3 and Latin 
ones over spatial coordinates. The metric signature is ( — 1-++) and we keep 
the discussion in four dimensions. 



2.1 Action and equations of motion 

We assume a gravitational action of the form 

- i a;(0) V^0V^0 -V {<!>) + ^ Xp, } , (5) 

where g is the determinant of the metric, k^ = SnG = STr/mp is the gravita- 
tional coupling, is a scalar field with potential V (so that it is a modulus of 
the theory when V = 0), fi and u are generic functions of 0, £gb is the GB 
invariant defined in Eq. (1), 04 is a constant, and £p. are the Lagrangians for 
fluids minimally coupled to 0. ^ 



Note that f2{4>) = — go'-^Cl^) iii the notation of [9]. 



The Einstein equations, SS/Sg^'^ = 0, read 
where 



f 1 

-4i?V^V,/2 - 8i?^,n/2 - 8i?(/^)V.V,/2 + 16i?,(^V'^V.)/2 , (7) 
and 



Here □ = VV^ = {— g)~^ ^"^ d^ {yf — g d^). In order to get the first equation, we 
have used the contracted Bianchi identities and the fact that the variation of 
the GB term in four dimensions is a total derivative (see, e.g., [20]). With a 
general function of the Riemann invariants in arbitrary dimensions one would 
get more (higher-derivative) terms. Taking the trace of Eq. (6) we obtain 

S = S^^ = cjV,0V"0 + 4V + (3n - i?) a + 4(i?n - 2R"^V^Vr)f2 
= T = T/. (9) 

The equation of motion for the scalar field is 



-4a4«:^[(V^0V^0)(V./2V'^0 + /2n0) + 2/2V^0V>V^V.0] = 0, 

(10) 

where differentiation with respect to the scalar field is indicated with a sub- 
script ^^. When 04 = 0, only the first line is left. 



2.2 Cosmological equations 



The Friedmann-Robertson-Walker line element with flat curvature in syn- 
chronous time is ds"^ = —dt^ + a?{t) dxidx\ where a{t) is the scale factor. The 
nonvanishing curvature tensor components read 



Ro^^o = a^iH^ + H), Rikki = -a^H\ i ^ k, (11) 

Roo = ~3{H^ + H), Ru = a\3H^ + H), (12) 

R = Q{2H^ + H), (13) 

where H = a/a is the Hubble parameter. When each fluid component is 
perfect, the energy- momentum tensor becomes 

T^u = Yl[{Pi+Pi)u^Uu+Pig^,u], (14) 

i 

where pi and Pi are the energy density and pressure of the fluids, and m'* is 
the unit timelike vector tangent to a fluid element's worldline. 

The Friedmann equation Sqo = Tqq is 

hH' + j\H + ^K'hH^ - a4«;V20' = y (ptot + ^ a;0' + V^) , (15) 

where Ptot = HiPii while Eq. (10) becomes 



o;(0 + 3if0) +V^ + ]- a;,^02 - 2Ah^^H\H^ + H) - ^ fiA'^H^ + H) 
+3a4«:^0'[4/2(0 + H<P) + /s,^^^] = 0. (16) 

We rescale time as 

t^T={t- to)Ho , 

where Hq = i/(to) is the value of the Hubble parameter today. Then the big 
bang is at Tj = —toHo, Tj < r < until today, and r > in the future. 
In the following, a subscript denotes quantities evaluated today, dots are 
derivatives with respect to r, and H is rescaled so that Hq = a'^da/dr = 1 
today. Defining 

Pco = SH^/k'^ , 

^i,0 = Pifll PcO ) 

i^ = K(j)/^/?> , 

U = V/p,,, 

and assuming that only matter and radiation contribute to the energy-momentum 
tensor, one has that Eqs. (15) and (16) become 



PhH'' + hH^ + hH-l a,Pf2v' = Pm + Pr + \oo^^ + U, (17) 
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+ |a4/90^[4/2(</3 + if0) + /2,^0^ 



where pm = VLmfi/o? and p^ = Qrfi/d'^- The former is a cubic equation in H, 
which can be rewritten as 



y{x) = bix^ + 62^^ + 63X - p = , (19) 



where 



x = H, (20a) 

bi=Pf2, h = fi, 63 = A, (20b) 

P = Pm + Pr + ^^^^ + U + l a^Pf2'P^ ■ (20c) 

Finally, the trace equation, Eq. (9), becomes 
u^''-W + h + ^Hh + 2{2H^ + H)h+l3[H^h + {'iH^ + '^H)Hh]=Pr^. (21) 

2.3 Theoretical models 



In the context of string theory, the field (/? can be interpreted as the massless 
dilaton field arising in the loop expansion of the low-energy effective action, 
with 



(22a) 
(22b) 

in the so-called 'string frame' and at tree level, where A = 1/4, 1/8, for the 
bosonic, heterotic and type II string, respectively. The coefficient 04 is fixed in 
order to recover the three-point scattering amplitude for the graviton [21,22]. 

Although we have set D = 4 from the very beginning, in general string and 
supergravity actions are defined in D = 10 (or 11) dimensions, and then com- 
pactified down to four dimensions (presently we ignore the braneworld case). 
Each compactification radius is associated with the vacuum expectation value 
of a modulus field. In the simplified case of a single modulus (one common 
characteristic length) and the heterotic string (A = 1/8) with stabilised dila- 
ton, we can identify (/? with such modulus and, following [1], we have 

/i = l, cu = 3/2, f/ = 0, (23a) 

/2 = -|eo ln[2e^V(ie^'')], «4 = 0, (23b) 

where 77 is the Dedekind function and .^0 is a constant proportional to the 
4D trace anomaly that depends on the number of chiral, vector and spin-3/2 



h = -^ = e-^^ 


f/ = 


/2 = ae-^^ 


a4 = -1, 



massless supermultiplets of the N = 2 sector of the theory. In general, it can be 
either positive or negative, but it is positive for theories in which not too many 
vector bosons are present. At large Iv^l, /2 ~ ^ocosh(v/3y9), an approximation 
which was shown to be very accurate even at small field values [23,24]. 

It is therefore reasonable to consider general functions 



/i = 7iiexp(7i2^), (24) 

/2 = 721 exp(722V?) + 723 exp(-722V5), (25) 

u; = 73iexp(732V?), (26) 

f/ = 74iexp(742V9), (27) 

and constrain the constants '-fij via the background and perturbed equations 
for real initial conditions. This will be done in Sections 4, 5 and 6. In particular, 
we shall see in Section 4.3 that the magnitude of the factors 7^2 in the dilaton 
exponentials is crucial to determine the quantum stability of the solutions, 
while for the modulus case it is not particularly important. 

The choice of potential, given by Eq. (27), is motivated by particle physics 
considerations. If we were to identify ip with the dilaton or modulus of a string 
model, then the most likely dynamics for these fields, leaving fluxes aside for 
the time being, would be nonperturbative effects. Those can come from either 
instantons or gaugino condensation [25] and are, in general, exponential in 
nature. 



3 Semi-analytic and numerical solutions 



In order to solve the Friedmann equation, Eq. (19), we introduce some useful 
quantities [26], 



x^ = -—- = ^, (29) 

3fei 3/3/2' 

If /i is constant, then r is positive and one can define h = \fr^ which is the 
distance (on the y axis) between the infiexion point and the local extremum. 
In general, for a cubic with real coefficients, one has only one real root when 



y'^ > r, three real roots (two or three coincident) when y1 = r,'^ and three 
distinct real roots when y1 < r. In the latter case h is well-defined. 

Let us write down r/y1 as an expansion in hi = /2/3 <^ 1, which is equivalent to 
an expansion in f3 if /2 is not fine tuned. The linear term vanishes identically, 
leaving 

;^ = 1 + ji (/' + 4p/i) {hpf + Oif3'). (31) 

Therefore the condition to have one single root is 

A' + 4p/i<0. (32) 

This is not possible if /i > 0, since p > for —a; < 0(1). If Eq. (32) holds, 
the root is Cardano's solution 

xi = x^ +P++P-, (33) 

where 

^^ = 2^ {^y* ^ "fy^) ■ (2^) 

In order to understand the behaviour of the solutions, we have to expand them 
up to 0{P). To lowest order in P, Eq. (33) becomes Xi ~ 3x* = 0{P~^), which 
is too large compared with today's value x ~ 0(1) and is not compatible with 
observations. This is true as long as /i, /2 ~ 0(1). 

The case f\ < is problematic in a Minkowski background, where the field ip 
acquires a constant vacuum value. In a FRW background, the reality of the 
solutions requires 
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l+4p^eR, (35) 



p2 

which, during radiation domination, becomes 

ifl^ - 2|/i|) ^' > 4|/i| [U + fi.,o(l + z)'] , (36) 

where z = a~^ — 1 is the redshift. This theory looks like a Brans-Dicke model 
with potential, for which there are quite strong cosmological bounds [27] at big 
bang nucleosynthesis (BBN). In other words, unless the potential is negative 
and fine tuned to be comparable to Qr^z^BN ~ 10^^, the field ip is subdominant 
to radiation. If so, Eq. (36) cannot be satisfied in general and at early times 
/i cannot be negative. This implies the existence of three real roots, of which 
we have to choose the positive one consistently with the initial data. As soon 
as the solution is chosen, the dynamics of the universe is set once and for all. 



We exclude this case a priori, since matter/radiation and the scalar field corre- 
spond to independent degrees of freedom. However, we are going to see that y^ ~ r 
approximately. 
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We conclude that a necessary (but not sufficient) condition to have a viable 
cosmological solution, compatible with BBN, is^ 

/i > 0. (37) 

The three distinct real roots of the cubic equation can be found via the La- 
grange resolvents method (see, for instance, [29]). One solution is given by 
Eq. (33), while the others are 



X2 = x^ + Cp+ + C'^P~, (38) 

X3 = x^ + C'^P+ + CP-, (39) 

where ( = (—1 + i\/3)/2 is a primitive third root of 1. Defining Cj = sgn(/j) 
and expanding to order 0{f3) we obtain 



X2 = -^-^{e,g + e^) + 0{f3), (40) 

xs=^-^{e,g-e^)+0{f3). (41) 

When €2 = +1, X2 is negative and x^ is positive, whereas when £2 = —1 the 
vice versa is true. In both cases, the positive solution is the same up to 0(/3^) 
and reads 



x^x^^^ +x^^^p, (42) 



x-' = ^ig-e,), (43) 



X 



(0) ^ 1/1 

2/1 _ 

(1)_ 4. 



2ffg 



(/i + p/i) 9 - if I + 3p/i) ei . (44) 



When 04 7^ 0, (7 must be expanded at the same level (then 04 = in x^^' ef- 
fectively). Note that x^^' is the solution of the quadratic Friedmann equation 
when f3 = 0. When /i = 1, Eq. (43) correctly reproduces the standard Fried- 
mann equation, x*^"-*^ = p. x^^^ is well defined in this case, x^^^ = — /2P/2. Note 
that there is no degeneracy in the roots, although they are formally given by 



^ In fact, one can find out Eq. (37) by promoting the number of e-foldings A^ = Ina 
to a scalar field, that is, a particle degree of freedom (sometimes this is referred 
to as Nordstrom theory [28]). Then, its kinetic term has a factor /i in front of it. 
Eventually one can find via the same trick that Q2 > 0, see Eq. (71). However, 
N interacts with ip and it is not clear in this approach how to define the energy- 
momentum tensor of each field. 
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the same expression once the sign of /2 is chosen. Therefore the cosmological 
problem is well posed and consistently resolved. 

Near a singularity, when the acceleration a is large enough, the approximated 
solutions for H diverge relative to the numerical ones, but we have checked for 
a number of examples that they agree at later times. In the cases described in 
Section 6, the agreement is very good for not too high redshifts. 



3.1 Numerical solutions and initial conditions 

Let us consider the Friedmann equation, Eq. (17). Taking its time derivative 
and using Eq. (18) we have 

And + Ai2V = Bi, (45) 

and, rearranging the equation for the scalar field, we obtain 

A21 d + A22^ = B2. (46) 

In the previous two equations the coefficients read 

Au = -{Pf2,^^H + f,), (47) 

a 

Au = f3f2,^H' + j\,^, (48) 

Bi = {3(3f2,^H^ + h,^)H^ - {Pf2,^^H^ + fi,^^ + Su) 0^ + 2f,H^ 

-(3p„ + 4p,)-|a4/3/20S (49) 

and 

A2i = --{Pf2,^H' + fi,^), (50) 

a 

A22=oo + la^Pf2^\ (51) 

B2 = fi,^H^ - SuHip - i u,^ v^2 -u^-l a^l3^\Af2H + /2,» . (52) 

Note that, in our units, both a and (/? are dimensionless. 

If the 2x2 matrix A with elements Aij is not singular (which is true in 
general, because the Friedmann equation and the equation for the field are 
independent), we can solve for a and (p. In fact, the determinant A = det A is 

A = An A22-A12 A21 = -[(2a;+9a4/3/20') {l3f2,^^H+h) + {l3f2,^H^+h^^f]. 

(Jj 

(53) 
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The equations can be written as X = A ^B, where X = (a, (/?)* and B = 
{Bi,B2y. Exphcitly, 

B1A22 — B2A12 , . 

a = , (54) 

y^= \ • (55) 

As we shall see below, the condition 

a > (56) 

at some time restricts the parameter space, so that not all sets of 7ij's (and, 
in particular, /2) will give rise to an acceleration era. The details of the model 
will depend on the choice for the reference time, which can be either today 
(late-time acceleration), near r, (inflation), or both. In the latter case there 
would be two constraints fixing the solutions, but then one would have to take 
into account the reheating phase in the early-time dynamics. ^ However, our 
choice for the time variable is not suitable for analyzing the inflationary phase, 
the interval between r, and Tbbn being too squeezed in our scale. Moreover, 
it seems difficult to reproduce the correct power spectrum in these theories, 
at least in the case where a collapsing phase precedes the big bang [30,31]. 
In general, the presence of a nonvanishing potential or cosmological constant 
helps in finding accelerating solutions. 

The set of initial conditions is Qrn,o ~ 0.25, fir,o ~ 8 x 10~^, a(0) = 1, d(0) = 1, 
and 0(0) = 00 as given by the solution of the Friedmann equation, Eq. (17), 
evaluated today, /3/2(v^o) + fiivo) + /i(v'o) = Po- 

|a4/3/20o + |^O0O-[/5/2,¥'(V^o) + /l,v>(V^o)]0O+[^m,O + ^r,0 + f^O " /l(V^o)] = 0, 

(57) 
which is quadratic in 00 for 04 = (i.e. the compactification modulus). In that 
case, when ujq > and the last coefficient is dominated by — /i < 0, there are 
typically two branches 0o^ < and (p^ > 0, roughly corresponding to €2 = +1 
(x3 root) and £2 = — 1 {x2 root), respectively; in general /2 can change sign 
during the evolution. In the dilaton case, cuq < and there are four solutions: 
if 04 = 0, (fQ is the larger dilaton root, which is positive. 

There is only one arbitrary initial condition left and one is free to choose it 
among our variables. The most attractive and physical choice might be today's 
acceleration, oq. However, in order to express (po as a function of cio one should 
invert Eq. (54), which can only be done numerically. Therefore we shall choose 
(p{0) = (/9o as the free initial condition. As one of other alternatives, one can fix 



^ Since the scalar field is assumed to be dynamical today, reheating should precede, 
for instance, via gravitational particle production. 
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(/9o and 00 and solve for Uq. We are also allowed to consider f3 in the interval 
< /9 < 1 without too strong a fine tuning. For an exponential /2, in fact, one 
can add a positive constant (/?* ~ O(IO^) to ip and define an effective coupling 
/3gg = /^e"^*. From now on and unless specified otherwise, we shall drop the 
subscript from PcS- 

The string nonperturbative regime occurs for y9 = 0. To impose ipQ <^ 1 
today would both be consistent with general relativity bounds in the case of 
the dilaton (the effective gravitational coupling being proportional to e''') and 
allow an evolution qualitatively different from GR in the case of the modulus 
with stabilized dilaton. However, we shall argue that this choice does not 
prevent the rising of instabilities. 

Since we want to see what the late-time behaviour is for these theories, the 
equations are integrated forward in time. It should be noted that a backward 
integration in time may not be faithful because it might lead to an attractor 
(which is actually a repeller if integrating from an early time up to now) that 
can be, in general, not consistent with a standard GR evolution. The numerical 
solutions are found to be well behaved and numerically stable, in the sense 
that for each run we plotted \(T,qq—Tqq)/(T,qq+Tqq)\ and |(S— T)/(S+T)|, and 
checked that the Einstein equations are satisfied by the numerical solutions at 
least up to one part to 10^. 

We conclude this section with a remark on the dilaton. Because of the negative 
sign of u, dilaton models are characterized by a contracting phase followed by 
an expanding one. At the minimum {H = 0), the scale factor must be smaller 
than at the nucleosynthesis, amin < Obbn ~ 10^^. This constraint drastically 
restricts the parameter space. The Friedmann equation at the beginning of 
the expansion is p = 0. Since Pm + Pr ^ lO'^^, one should impose, neglecting 
the term in a4/9, that |c<j0^/2 + U\ ~ \ujip'^\ > 10^^ at aminj a condition hard to 
achieve. Anyway a > in this class of models, and we can rule it out. 



4 Cosmological perturbations and no-ghost constraints 



One of the key ingredients of this work is the non-minimal coupling between 
the scalar field and the gravitational sector. Such interaction has been intro- 
duced in many different contexts and, probably, scalar-tensor theories repre- 
sent the simplest example. String models typically couple the modulus field 
not only with the Ricci scalar, as usual in a scalar-tensor theory, but also with 
the Gauss-Bonnet term: the hope in allowing for a runaway field is that it 
may also explain today's acceleration of the universe. 

It is well known that, for standard scalar-tensor theories with a coupling of 
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the form i?, the scalar field itself cannot acquire negative expectation values, 
otherwise the graviton would become a ghost. By expanding the action at 
second order in the perturbations on a Minkowski background, the kinetic 
term of the graviton tensor modes is 0o h^j ^h^i for a constant 0o- Therefore 
it is required that 0o > 0. In the case of a coupling with the GB term, the 
expansion about a Minkowski background is trivial in four dimensions. This 
is because 0jCgb becomes a total derivative, being = 0o in the Minkowski 
vacuum, and the GB term automatically disappears from the equations of 
motion. 

The same effect is obtained in any other background if the dilaton is stabilised 
by some mechanism. However, in the presence of a time- dependent dilaton 
things become more complicated. It should be emphasized that this is exactly 
the scenario discussed in this paper, where the background is not Minkowski 
but Friedmann-Robertson- Walker. Typically, in these backgrounds the field 
acquires nontrivial dynamics and one expects that a coupling with the Gauss- 
Bonnet term may lead to inconsistencies of the theory, such as a change in 
the sign of the kinetic terms of some physical degrees of freedom. To see more 
in detail what we mean by this, we consider the gravitational perturbations 
about a FRW background, using the standard technique of studying the tensor, 
vector, and scalar contributions separately [32,33]. 

It is convenient to work in conformal time f] = J dt/a. The linearly perturbed 
fiat, FRW line element can be decomposed as follows 



ds"^ = a^{r]) {-(1 + 2$) dr]'^ + (2 5 ^ - S,i) drf dx' + [hy + 2Fij 

+ (1 - 2^) 5ij + 2E^ij] dx' dx^} , (58) 

where hij is a 2-tensor with spatial indices. Si and Fi are vectors, and $, ^E', 
B, and E are scalars. 



4.1 Tensor perturbations 



The FRW line element at zero curvature is conformally fiat and it is possible 
to choose the transverse-traceless harmonic gauge for the tensor perturba- 
tions [34] 

h^o = = V. djh/ = 0. (59) 

Ignoring the other modes, it is long but not difiicult to find that the action at 
second order in hy is 
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5(2)5tt= I dr]d^xa^ 



h^h^'-'i-""'^^ 



hh^h'j 



+ \{f!^-2nfl)h^.h'n + ..., (60) 



a2 



where □ = —d'^ + d'^ is the Minkowskian d'Alembertian, primes denote deriva- 
tives with respect to r], Ti = a' /a, and we have kept only terms containing 
second derivatives of the perturbation h^j (the dots stand for hnear-derivative 
terms, no effective potential). The integrand can be written as 

-[K^{r]) - K^iv)] hhdlh', + KM h^^^lh', , (61) 

where 



^i(^) = ^ + ^(/2-^/2), (62) 

K2{v) = \{f^-2nf^). (63) 

By performing a change of the time variable t] ^ fj, we have that 

7-\ 2 



^^^ ^2, , ;, , d^V 



Therefore if 

p-1 , (65) 



dv I K^ _ n-l 



dri \ Ki- K2 
the second-derivative operator becomes 

I df] d^x a^ P{f]) Ki{f]) hU ah'j + ... . (66) 

We require the theory to be ghost free. Combining (65) and (66), one has a 
ghost if -ft'i < and, because of Eq. (65), for Ki — K2 < 0. In other words the 
kinetic term of the field is well defined only if 

Ki>0 and Ki - K2 > . (67) 

In universal time the above conditions can be written as 

/i + 8«:2/2>0, (68) 

h + 8K^Hf2>0, (69) 

or, in terms of r and rescaled units, 

TTl: Q^ = f^+f3f^>0, (70) 

TT2: Q2 = fi+PHf2>0. (71) 
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Incidentally, the speed of propagation Stt of tensor perturbations is indeed 
equal to QxjQ-i- This speed needs to be positive definite, otherwise the system 
becomes unstable: the second-order differential operator in Eq. (66) would 
become (plus or minus, depending on sgn(i^i)) a 4D Laplacian rather than 
a d'Alembertian and perturbations would grow exponentially. Therefore this 
kind of instability occurs if either TTl or TT2 are violated, whereas a ghost 
appears when both these conditions do not hold at the same time. 



It is important to notice that the conditions TTl and TT2 still hold in the 
presence of more than one scalar field, because the scalar degrees of freedom 
do not mix with the tensor ones: for example, j^ may be a function of, say, two 
fields which can be thought of as representing the real and imaginary parts 
of a complex scalar (see Section 7). These conditions should hold at all times 
and a numerical simulation is necessary to make sure that both Eqs. (70) and 
(71) are satisfied during the evolution of the universe. It is also clear that, if /9 
is small and /2 and its derivatives are always < 0(1), then the equations are 
satisfied providing that /i > 0, which is the no-ghost condition for a scalar- 
tensor theory deduced already from the background equations, Eq. (37). Of 
course, for a de Sitter background with constant (/? both conditions hold but, 
in general, this is not the case (in particular, a dS background is actually 
achieved also when (/? oc r). 

We can already restrict the parameter space for the initial conditions by requir- 
ing that both TTl and TT2 should be valid. As an example, for the dilaton, 
TT2 implies 



l-^/3v3o>0. (72) 



For values of /9 ~ 10~^, TT2 is not satisfied, giving rise to an instability 
when, for example, (/9o ~ 10 and we choose the branch >^q ~ O(IO^) (f/ = is 
assumed). The level of fine tuning in this example is still modest, but it is clear 
that it corresponds to an unphysical situation, since the effective gravitational 
coupling would be strongly enhanced today. 

It is interesting to note that the dilaton case can in principle have ghost 
modes on FRW, although this requires a tuning of the parameters in contrast 
with experiments. In fact, in these theories, ^o is typically related with the 
time derivative of the Newton constant, that is |7i20o| ~ \B.q^ G^^ dG j dt\Q < 
1.75 X 10~^ (see [35]). Then it is true that viable string-inspired cosmologies, if 
any, are ghost-free like the mother theory, but this property is not automatically 
satisfied by construction. This result is independent of the value of a^. 
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4-2 Vector perturbations 



It can be proved that, in the absence of an anisotropic fluid and introducing 
the quantity Vi = Si + aFi, the following quantity is conserved 



a^fi + 8 k'^ H f2) V(^ij) = const. 



(73) 



This is equivalent to saying that vector perturbations do not propagate like 
waves or, in other words, the angular momentum of the fluid is conserved 
(see [8,10]). The other equation for the vector perturbation in vacuum states 
that Vi is an harmonic function satisfying the Laplace equation, djd^Vi = 0. 



4-3 Scalar perturbations 



If we can neglect the matter contribution, i.e. at late times or in vacuum, the 
calculations simplify considerably by choosing the uniform field gauge 6(f) = 0. 
In this case it is possible to show that the action for the potential ^f can be 
written as [81 



^^'^'5sc = 2 I d'xdta:'Qsc 



^2_ 



'SC 



{d,^f 



where 



and 



Q 



SC 



3 (fi+SK^ H^ f2)^ 
2k2 f^+8K^Hf2 






^ 2 fi+SK^HfiJ 



(74) 



(75) 



ssc = 1-4 









H'^HM^§Jfy2a,n^h, 



, , i2 I 3 (/i+8k2H2/2)2 ^4 f i4 



(76) 



Therefore, it is clear what the no-ghost conditions for the scalar modes are. 



Qsc > and Sgc > 
The first condition can be rewritten as 

3{j\ + 8K'H'f2'' 



or 



where 



K^UJ 6^ + 



SCI : 



2 f\ + 8K^Hh 



+ 12a4/€V20 > 



q = uoLp + 



2Q 



27 
' + — a4/3/20'>O, 



Qs = fi + PH'f2 = H{Q2-h) + h. 



(77) 

(78) 

(79) 
(80) 





If cj > and a^ = 0, this condition does not add anything new to the tensor 
constraints. 

The second condition, that is trivially satisfied if /2 = and a^ = 0, is 



SC2: s,r = l-^ h\^] \^-H- 

>0. (81) 

This is a new constraint, being independent from the others we have already 
encountered. Therefore there are, in general, four different no-ghost conditions 
that should be satisfied in order for the theory to be consistent and free of 
instabilities. ^ Note that s^c is interpreted as the speed of propagation of the 
perturbation. 

The no-ghost conditions do not depend on the scalar potential. This is because, 
intuitively, the sign of U determines whether (/? is a tachyon (which is an 
instability quite recurrent in FRW), but does not affect the kinetic term. In 
any case, as we shall see below, the potential does modify the background 
solution and, indirectly through the no-ghost constraints, their stability. 

When /2 is constant or (3 is very small, Eq. (79) gives a bound on the value of 
the coefficients 7i,- which does not depend on ip. Assuming 732 = 712 and 711 
as in the dilaton case, one can factorize the exponentials and get (/? = 0) 

2— + 7i'2>0, (82) 

7ii 

where we have imposed the TT condition 711 > 0. This equation is valid in 
the GR limit, which must be reached at some point during the evolution of 
the universe. For this reason, actually, it holds at all times as a necessary (but 
not sufficient) condition. For a minimally coupled scalar field (712 = 0) and 
positive gravitational coupling, this condition is trivially 731 > 0, while for a 
dilaton-like field 711 = —731 = 1, it implies that I712I > V2. The string dilaton 
respects this bound. 

Defining 

0BD = TT e'^"'^ , c^BD = ^ , (83) 

2 7i2 711 

the Lagrangian with f3 = reduces to a Brans-Dicke theory in Einstein gravity. 



^ Strictly speaking, not all violations of these conditions will lead to ghost instabil- 
ities, as explained below Eq. (71). For convenience, we shall keep calling the above 
conditions 'no-ghost'. 
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£bD = 0BDi? - ^V^^BoV^^BO - ^(0Bd), (84) 

0BD 

and Eq. (82) becomes the usual no-ghost constraint cUbd > —3/2 (together 
with the TT condition 0bd > 0) for Brans-Dicke theories. 

Note that all the results in the scalar sector are obtained in a gauge-invariant 
framework [10]. ® 



4-4 On superluminal propagation of perturbations 



There is another issue, related to the speed of propagation of the perturbations, 
which can further constrain cosmological solutions. Particles that propagate 
at speed faster than light can generate paradoxical situations in which Lorentz 
invariance is broken and preferential frames in which the physics is well de- 
fined are undemocratically selected. In Minkowski background, superluminal 
propagation is associated with breaking of causality, ill-posed Cauchy problem 
[36] , the possibility of closed time-like curves, and violation of the null energy 
condition. The latter is usually responsible for instabilities, although there are 
stable systems in which the null energy condition does not hold [37]. Superlu- 
minal particles are also associated to the spoiling of conventional black hole 
dynamics and thermodynamics [38,39]. 

To extend this analysis to curved backgrounds is not easy and, to the best 
of our knowledge, only partial evidence in support of bad causal behaviour 
has been provided so far [40,41,38]. For instance, the presence of superlumi- 
nal modes would imply that for any comoving observer on a inertial reference 



Actually, in the more general case where the term /i((^)i? in the action is replaced 
by f{4'iR) &!■ some smooth function /, one can find two scalar degrees of freedom 
obeying independent Mukhanov equations. This can be understood by recalling that 
f{R) gravity has a hidden scalar degree of freedom, which adds to that coming from 
the GB term. However, one can see that there is always a unique set of no-ghost 
conditions as follows. 

The two scalar modes are governed by Eqs. (104) and (105) of [10] (see the orig- 
inal paper for the notation; here C2 = C3 = 0), and the former can be writ- 
ten as {a^Qsc)~^d{a^Qsc^)/dt = Ssc(5j9*<l>)/a^. Equation (105) can be recast as 

{a^Q^)-^d{a^Q^,^)/dt = Ssc{did'^)/a'^, where § = '^aQ2/[H + Q3/{2Q2)] and 
Q^ = a^'^{QscSsc)~^- After requiring Ssc > 0, to impose either Q^ > or Qsc > 
results in precisely the same condition. 

In the case of study (no anisotropic stress, Einstein-Hilbert leading term) there 
is only one independent scalar field which can be canonically quantized, provided 
the speed of propagation is real. Consistently, in the uniform field gauge {5(p = 0) 
$ = *. 
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frame (say, an experiment to detect gravitational waves), the passage of a su- 
perluminal signal would be felt as a breakdown of causality. On the other hand, 
the Mukhanov equations for the perturbations are defined on a conformally 
fiat background, and this alone might be sufficient to justify the imposition of 
the suh-luminal (SL) constraints 



SLl: s^^ = Qi/Q2<l ^ f2-Hf2<0, (85) 

SL2 : Ssc < 1. (86) 

Our perspective is to dismiss as unviable solutions on which perturbations 
propagate faster than light. 

Note that, when 04 = 0, Eq. (81) can be rewritten as 

6q%{l - s,c) = -(1 - s,,) +4e(l-^], (87) 



QV ' ' ' \ Q 



3/ 

where e = —H/H^. In a non-superaccelerating universe (e > 0) in which the 
SLl and TT2 conditions are satisfied, the first term is nonpositive, while the 
second one is non-negative if, and only if, 

Q3>/i ^ /2>0. (88) 

Equation (88) is necessary to guarantee the sub-luminal condition SL2, 1 — 
Ssc > 0. 



4-5 When other fluids are present 



All the above results were found in the absence of other perfect fiuids in the 
action, and it is natural to ask what happens in more realistic situations when 
both pressureless matter and radiation evolve together with 0. 

At late times both matter and radiation are subdominant relative to (re- 
garded as dark energy) and they can be neglected in an asymptotic solution. 
Therefore, if a solution does not satisfy the ghost constraints at late times, it 
is sufficient for us to discard it as unviable altogether. 

Problems may arise when trying to check a solution which is ghost-free asymp- 
totically in the future, or when one assumes that dark energy is not the mod- 
ulus but another fiuid with w ^ —\ today. However, we do not expect the 
four ghost conditions to be modified in the presence of extra perfect fiuids, 
provided these are minimally coupled to both gravity and the scalar field. 
The tensor constraints would be unchanged since there are no tensor modes 
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in isotropic fluid perturbations. On the other hand, to any new fluid compo- 
nent would correspond another scalar mode to disentangle from that of the 
modulus perturbation (for a multi-field example, see [42]). 

There is evidence that the constraint equations already written would not be 
affected. A way to see this is to note that the usual no-ghost constraint cJqd > 
—3/2, Eq. (78) for Brans-Dicke theories in Einstein gravity, comes also from 
a conformal transformation from the Jordan to the Einstein frame, where the 
(no-)ghost mode becomes apparent. This constraint does not depend on the 
frame choice, nor upon the presence of non-minimally coupled perfect fluids. 
In the moduli case, if Eq. (78) were modified when a new fluid component is 
taken into account, the only way to recover the Brans-Dicke case would be to 
have couplings like 6(/2, /2) c{p,p) for some functions b and c. But at the level 
of the background action such couplings are absent, and in the ghost conditions 
only couplings already present from the beginning can appear (for instance, 
both /i and /2 couple to gravity, while u does not). A direct inspection of 
the perturbed Einstein equations is difficult for assessing this result more 
rigorously and, for the moment, we shall regard it as a conjecture yet to be 
proved. 



5 Ghost conditions for models in literature 



Although we have considered only inhomogeneous (i.e. cosmological) pertur- 
bations, it is common in literature to study the attractor solutions of modified 
gravity models via a phase-space analysis. A natural question one may ask 
is: can an attractor, stable against homogeneous perturbations, be unstable 
against cosmological perturbations? The answer is yes, and we shall provide 
several examples. The inequivalence of the two approaches was stressed also in 
[19,43]. A weaker (and maybe more obvious) result holds, namely, phase-space 
unstable solutions can be free of ghosts. 

We can compare the prediction of the no-ghost constraints with some known 
cosmological expanding solutions [14], with Gauss-Bonnet parametrization 
in four dimensions. To make contact with past notation, we switch back to 
synchronous time and set /t^ = 1. 



5.1 Asymptotic solutions and modulus scenarios 



At late times one can identify three different limits: a 'low curvature' regime 
in which the Einstein-Hilbert action dominates over the Gauss-Bonnet term, 
a 'high curvature' regime where the quadratic term dominates over the lin- 
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ear one, and an 'exact' regime in which both terms contribute at the same 
level. Asymptotically, one can consider a power-law expansion, a ~ t'^^ with 
logarithmic modulus (/i = 0) 

i^~-, 0=-, /2~^t-, (89) 

t t 8 ^ ^ 

where ctj and .^o are determined by which of the three limits above one imposes 
to reach. Any fluid content has been damped away at earlier times. ^ For the 
compactification modulus (/i = 1, a; > 0), the no-ghost constraints read 



TTl : 1 + ^oa2{cr2 - 1) t"^"' > 0, (90a) 

TT2 : 1 + ^0^1 f^2 1"''-^ > 0, (90b) 

gC2; a|^ (eo^,V2t--3)2 (eoa2)M 



t2 1 + eofTias t-2-2 3(1 + ^^aia2t''^~^y 

4^2(<72-3) 
0-2-0-1 



X 



+ 3)t3'"2-8 



60-10-2 



> 0. (90c) 



Here we do not discuss the viability of these solutions as models of dark energy, 
and limit the comparison to the stability issue. After looking at attractor 
solutions, we shall see a couple of numerical examples. In this section only, a 
solution is said to be stable or unstable in the phase-space sense. 



The low-curvature solutions are o"i = 1/3, (72 = ±2/3, and both branches are 
stable. ^° The TT and SC conditions are trivial: as regards SC2, the leading 



term is a^t ^ > 



There is also a high-curvature asymptotic solution for the modulus with zero 
potential, which requires a nonvanishing cosmological constant [14]: o"i = 1, 
a2 = 4:, V = 12^. Such solution is an attractor only when ^ > 0. Both TT 
conditions require ^ > 0. However, SC2 gives 



'SC 



-46 >0 ^ 6<0. (91) 



Therefore the attractor solution has always ghosts. Conversely, non-attractor 
solutions (6 < 0) are always ghost-free. According to our interpretation of the 
ghost conditions, one can infer that it is never possible to achieve a physical 
high-curvature regime asymptotically. 

Figure 1 shows the behaviour of e^^ and (p{eH)~^ and the convergence to 
the attractor at e^^ ~ Ht ~ ai and (p{eH)~^ ~ 0t ~ o"2- The condition 
SC2 is not satisfied at late times. The numerical analysis unravels a hidden 



^ In the notation of [14], Uj = tOi, ^0 = 45. Also, we do not consider solutions with 
sudden future singularities. 
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If one allows for an extra perfect fluid, it must be stiff matter, p ■ 
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100 



100 120 



Fig. 1. The modulus high-curvature attractor, for /3 = 10 '^, (/)o = 6.8, and (po ~ 0.7. 
Top: e~^ ~ o"i; centre: <p/{eH) ~ o"2; bottom: the ghost constraint SC2. 



issue for the asymptotic solutions. The ghost constraints Eqs. (90) are vahd 
only at late time and do not say anything about the early-time behaviour 
of solutions approaching the power-law attractor. Then one might encounter 
unstable modes even before the attractor is reached. This means that, of all the 
possible real initial conditions, some are excluded by the presence of ghosts, 
even those that would have led to an attractor in the phase space. 

Finally, an exact unstable solution is cxi ~ 0.21, a2 = 2, ,^o = 54.12 [1]. The 
TT conditions hold but SC2 does not, the latter being Sgc ~ —2.29 < 0, and 
the solution is unstable also against inhomogeneous perturbations. 
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5.2 The Nojiri-OdintsovSasaki modulus 



This model for dark energy ([13], hereafter NOS) is defined by making the 
choice 



/i = 1, cu = ±1 , 



(92) 
(93) 
(94) 



where /o, a, 0i, Vq are all constants. By choosing an appropriate value for Vq it 
is possible to have an accelerating universe today. Let us study the presence of 
ghosts in such a model. It is useful, before performing a full numerical simu- 
lation, to have a look at the initial conditions which determine the parameter 
space for which there are ghosts. In other words, the conditions TTl, TT2, 
SCI, and SC2 can be studied for today's initial conditions, which are functions 
of just 00- 

For a choice of the parameters of order unity {P = a = (pi = fo = 1, u = 1, 
and Vq = 0.7), the reality of 0o requires that 0o ^ —0.3 for both the Re- 
branches. In this case it is possible to check the no-ghost conditions (again, 
SCI is trivial) and we find that they are all satisfied but SC2 (negative squared 
propagation speed, see Fig. 2) and SLl. This is consistent with the numerical 





0.2 0.4 0.6 0.8 1 

00 

Fig. 2. No-ghost condition SC2 for initial conditions 0o in the NOS model. Left: 
branch; right: (pt branch. 



analysis made by choosing the negative branch and 0o = 0.2. In this case there 
is a late-time power-law attractor with ~ Int, and both SC2 and SLl are 
not satisfied (see Fig. 3). 



25 




100 







2.5 

2 

1.5 

StT) Ssc 1 


\ 

V : 




0.5 



-0.5 


/ 


J 



0.001 0.01 



0.1 1 

T 



10 



100 



Fig. 3. Late-time behaviour for e ^ ~ cJi (top), (j)/{eH) ~ o"2 (central), and the 
speed of propagation for the scalar (solid) and tensor (dashed) modes (bottom) in 
the NOS model. 

5.3 Asymptotic dilaton scenarios 



The logarithmic dilaton with no potential [44] has a stable, low-curvature 
asymptotic solution ai = l/VS ~ 0.58, (T2 = a/3 — 1 ~ 0.73 (there is no 
expanding solution in the presence of a cosmological constant). Terms in /2 can 
be neglected and it is easy to see that the no-ghost and subluminal conditions 
are trivially satisfied. 
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5-4 Linear dilaton 



An exact de Sitter solution of the equations of motion can be obtained when 
(f) = vt and V = 0, a4 = —1, with v ^ 1.40 and H ^ 0.62 in the bosonic 
case (there is no such real solution for 04 = 0). For the GB parametriza- 
tion, this solution is stable. The no-ghost tensor conditions are satisfied, TTl 
being trivial and (^2 oc 1 — Hv ^ 0.13 > 0. SCI is satisfied, too, but SC2 
is violated, being s^c ~ —7.61. In [9] the instability arising when s^c < 
was already noted, and interpreted as a breakdown of the linear theory. The 
same authors claimed also (only as a preliminary result) that the inclusion 
of higher- derivative terms in the action would flip the sign of Ssc- However, 
the linear expansion of the metric is typically taken as sufficient to determine 
the viability of a noneffective quantum field theory from the point of view of 
unitarity, which is also our perspective. 



5.5 Cyclic scenarios 



Another case in which the action (5) plays an important cosmological role is in 
cyclic or bouncing scenarios, in which a contracting phase before infiation can 
leave an imprint into the cosmic microwave background (CMB) perturbations 
(e.g., [30] and references therein). 

As an example, we take a dilatonic scalar field with 04 = and a potential 
corresponding, in the Einstein frame, to a negative cosmological constant, 
A = 1/8, 741 = 7j2 = —1- The scale factor and the condition SCI are shown in 
Fig. 4. We have checked that |(S-T)/(E+T)| < 2x 10"^ and |i//(Soo+Too)| < 
10^^^. Each bounce is nonsingular in H, and there is a correspondent cyclic 
increase of |g|. 

We notice that linear theory may break near each bounce [45], and the ghost 
constraints above might not be too well motivated in this case. However, a is 
never singular and g < always, even far from the bounce. The presence of 
a spin-0 ghost in the general class of bouncing models was already noticed in 
[31], where such instability is regarded as the key ingredient for the bounc- 
ing mechanism (however, that setup is rather different from ours). Also, this 
particular example is expected to be unstable because of both the negative 
sign of the potential and the choice for 7j2 (although later on we shall see 
that instabilities appear anyway for most choices of these coefficients). Here 
we show it just as another situation in which a non-minimally coupled scalar 
with a tailored potential and arbitrary initial conditions is plagued by instabil- 
ities. We checked that the background solution does not change qualitatively 
when 04 = —1 (the root closest to ipo was taken and found to be almost the 
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Fig. 4. The scale factor and the function q of the dilaton cychc solution described 
in the text, for (3 = 10"'^, (po = 1.8, and 0o ^ 1-2. The ghost constraint SCI, which 
is shown in one cycle, is never satisfied. 

same). The ghost constraint SCI is worsened by the extra term in 04, which 
is negative definite. 



5. 6 Homogeneous vs cosmological perturbations 



To summarise, we have given examples of stable phase-space solutions which 
have ghosts or other tree-level instabilities. Looking at the linearly perturbed 
equations of motion in phase space [14], one can see that the characteristic 
equation for the mass matrix has a structure which cannot reproduce that 
of the no-ghost constraints. While a phase-space analysis can track solutions 
which are stable in the classical sense, the structure of cosmological perturba- 
tions is in closer relation with stability at the quantum level, and the two are 
independent in all respects. 

So far we have limited the discussion to an action with the GB parametriza- 
tion, but solutions (attractor or not) are known also in the general case where 
the Riemann invariants R^, P, and Q have arbitrary coefficients [14]. The 
perturbed equations of motion would be very complicated but it is clear that 
there would be always a spin-2 ghost in the spectrum of the theory (this is 
trivial in a Minkowski background). In fact, for the tensor modes in FRW, the 
Lagrangian 

/i(0)i? + f2{(p){a^R^ + a^P + a^Q), (95) 

would give rise to a combination /iD/ij-' + ^(aj, /«, /«, H, . . . )D'^hi^ for some 
function u, which has always a ghost unless either /i or u vanishes identically 
at all times. ^^ While /i 7^ is phenomenologically necessary, the function u 



^^ Higher-derivative terms are usually responsible for violation of unitarity. Taking 
the example of a scalar field with equation (□ + 'yn'^)(p = and calculating the 
propagator G{—k'^) in momentum space, one sees that —G{—k'^) = k~'^ — (/c^ — 
7~^)~^. Then the particle content of this model is a massless scalar and a ghost 
scalar with mass 7~^, the relative — sign between the propagators being the origin 
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cannot vanish for arbitrary choices of a^ and /, (when /j = 0, m oc 02 + 4a3 
and the coefficients a, are fixed to the GB choice [7]). Then all cosmo logical 
solutions of the above generalized quadratic Lagrangian, Eq. (95), have ghosts, 
even if there are attractors. 

One can conclude that the studies of the homogeneous and cosmological per- 
turbations are inequivalent. 



6 Towards observationally and theoretically viable solutions 



So far we have shown that the presence of instabilities depends on the bound- 
ary conditions (hereafter 'b.c.') selected for each model. It is now time to ask 
whether there exists a nontrivial region in the parameter space of the b.c. 
such that the resulting evolution is free from instabilities and compatible with 
supernovae observations. The task is hard since the no-ghost and sub-luminal 
constraints are time dependent and they should be evolved for every b.c. and 
checked to hold at any time throughout the whole evolution. Here we do not 
fully address this problem but note that a necessary although not sufficient 
condition to have experimentally viable and (quantum mechanically) stable 
solutions is to impose the no-ghost and sub-luminal constraints, together with 
the present bound for a, to be satisfied today. Of course this does not guar- 
antee ghost freedom always, but we will see that the parameter space shrinks 
considerably nonetheless. 

In our particular class of models the parameter space / consists of 17 elements, 
I = {P, 'jij, 04, ipo, (pQ (/3o, Oo, Co, Co}. Of these, two are fixed by default (cq = 
1 = do in our units), while do, v3o, and tpQ are given by Eqs. (54), (55) (evaluated 
at r = 0) and (57), respectively. We can then specialize to the modulus case 
with arbitrary potential, leaving the effective GB coupling unspecified. Hence 

{711, 712, 721, 722, 723, 731, 732, 04} = {1, 0, ±1, VS, ±1, 3/2, 0, 0} . 

There is still freedom in the choice of one of the two branches for 0o . We shall 
take the positive branch (pQ when /2 > (and vice versa), the negative one 
being unviable because it is then much harder to satisfy the TT constraints 
on sufficiently long (/9o-intervals. 

The parameter space is then reduced, for each branch, to a two-dimensional 
submanifold, I = {(3,ipo}. The other parameters are either fixed by hand or 
extracted from the equations of motion. 



of negative-norm (or negative-energy) states. See, e.g., [46] for a short introduction 
to ghosts in higher-derivative theories. 
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In order to show the impact of each set of equations, we shall adopt a gradual 
approach by first imposing only the four no-ghost constraints (TTl, TT2, SCI, 
SC2), then the no-ghost and subluminal constraints (SLl, SL2), and finally the 
no-ghost and sub-luminal constraints together with the observational bounds 
for ao [47] 



2a: 0.22 <ao< 1.11, (96) 

la: 0.55 <ao< 0.97, (97) 

for a total of maximum 7 constraints. ^^ In principle, one may take into account 
even constraints on 'a, as in Table 5 and Fig. 5 of [47]. Here we will not need 
to consider this extra bound anyway. 



6.1 Modulus case 



The simplest available example is a modulus in strict sense, that is, with fiat 
potential, [/ = 0. In Fig. 5 all the constraints are shown for a particular value 
of 13 (see the caption for details). One can then project the intervals, along 
the ifQ axis, into which some or all the depicted functions are non-negative, 
and repeat the process by varying (3. The region in the ipQ-\nf3 plane where 
only the no-ghost conditions hold is shown in black in Fig. 6 (left panel). 
By imposing also the SL conditions, the allowed region is appreciably reduced 
(right panel), with an upper bound at /3 ~ 0.30. This reduction comes actually 
from SLl only (tensor modes do not propagate faster than light), being Ssc < 1 
at all times. However, the universe does not accelerate today, oq < 0, and the 
model can be discarded. To get acceleration, we must introduce a nontrivial 
potential. We can take an inverse exponential, U = e~^^^, as in the NOS case 
(see Fig. 7). The positive branch ip^ is real when ipo > 0.167 or ipQ < —7.976, 
but in the latter case it is easy to see, with a similar analysis, that the no- 
ghost constraints cannot be satisfied at the same time and there are always 
instabilities. 

Again, SL2 is trivial and the allowed region in the parameter space is mainly 
reduced by SLl (Fig. 8), with P < 0.34. When including Eq. (96), only a very 
narrow strip survives the selection, whose thickness depends on which sigma- 
level bound one assumes (Fig. 9). Note that the selection is mainly carried out 
by the lower bound for oq. In fact, /J > 10^^ is in the allowed region up to a 
maximum acceleration cq ~ 0.41, so that the model is rejected at the la level. 



-"^^ The bounds Eqs. (96) and (97) come from the joint analysis of SNLS [48], 
SDSS [49], WMAP3 [50], and 2dF [51,52] data, assuming a nonconstant effective 
barotropic index for the dark energy equation of state. 
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Fig. 5. No-ghost, sub- luminal, and acceleration conditions for a modulus field with 
[7 = and (3 = 10^^ as functions of the initial condition ip^ (positive branch). Solid 
curves with increasing thickness are Qi, Q2, q, and s^c- The dashed, dot-dashed, 
and dotted curves are 1 — Stt, 1 — Ssc, and 'do, respectively. 



\nf3 





Fig. 6. Region (black) in the (/Jq-Ih /3 plane satisfying the no-ghost (left) and no-ghost 
-|- sub- luminal conditions (right), for the positive branch if^. Here, U = 0. The 
model of the previous figure is at ln/3 ~ —6.9. 

Then, when Sq is increased by 0.01, there is a leap of 11 orders of magnitude 
in the GB coupling and the allowed region lies below f3 ~ 10"^^. 

We can learn a few things from these pictures: 

• As one would expect, when P decreases and the theory goes towards the 
GR limit, the allowed intervals on the (fQ axis for fixed (3 expand. 

• An interesting result is the presence of upper bounds for j3, above which 
there is at least one moment during the evolution of the universe when 
instabilities appear. This might mean that the GB Lagrangian, regarded as 
a perturbative expansion in [3, breaks down and higher-order corrections 
are necessary for consistency. 
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Fig. 7. No-ghost, sub- luminal, and acceleration conditions for a modulus field with 
U = e^^^^ and /? = 10"'^ as functions of the initial condition ip^ (positive 
branch). Solid curves with increasing thickness are Qi, Q2, Q, and Sgc- The dashed, 
dot-dashed, and dotted curves are 1 — Stt, 1 — Ssc, and 'do, respectively. 



ln/3 _ 




Fig. 8. Region (black) in the v^o-ln P plane satisfying the no-ghost (left) and no-ghost 
+ sub-luminal conditions (right), for the positive branch (pQ. Here, U = e~^'^'^. The 
model of the previous figure is at ln/3 ~ —6.9. 



In at least one example the sub-luminal condition, SLl, removes most of the 
parameter space. From Fig. 5 it is clear that, keeping only the solid lines, 
the allowed interval along ipQ is rather larger than when one imposes also 
the SL constraints. This is confirmed in the left plot of Fig. 6, where the 
positive branch allows also negative values for ip^. Here there is no upper 
bound for f3, which indeed means that it comes from the requirement of 
sub-luminal propagation. 

The condition of Oq is by far the most stringent bound for the parameter 
space. In general, not all choices for 7jj are compatible with observations. 
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0.24 0.26 0.28 0.3 0.32 

^0 



Fig. 9. Region (black) in the c^o-lii/3 plane satisfying the no-ghost, sub-luniinal, 
and acceleration conditions for the positive branch (p^ and modulus potential 
[/ = e"^^. The 2a level of [47] has been used, Eq. (96). 




Fig. 10. No-ghost, sub-luminal, and acceleration conditions for a modulus field with 
/2 < 0, {/ = e"^^*^ and (3 = 10^^ as functions of the initial condition (pQ (positive 
branch). Solid curves with increasing thickness are Qi, Q2, q, and Sgc- The dashed, 
dot-dashed, and dotted curves are 1 — Stt; 1 — Ssc) and oq, respectively. 

The above features are qualitatively the same when /2 < 0. ^^ In this case, the 
negative branch (^q is the viable one; the example U = e~^^'^ with the sign of 
/2 flipped is shown in Fig. 10. Again, the allowed region is for P < 0.30. 

Let us take ipo = 0.25 in the example of Fig. 7 (so that Hq fs 0.32). At late 
times one recovers an asymptotic evolution as given by Eq. (89) with o"i ~ 3 
and (72 ~ 2/v^ ^ 1.15 (see Fig. 11). With (T2 = 2/v^ the contribution of 
the potential goes as the scalar kinetic energy and the Hubble parameter. 



^^ We have checked that they do not change also for other values of 7j2, although 
the allowed parameter space changes and may be constituted, for each /3, by the 
union of disconnected intervals along (/9o-slices. 



33 







1.2 


3 




1.1 
1 




^-^ 


12.5 


/ 


^ 0.9 


2 


/ 


0.7 
0.6 


1.5 


/ 





3 

r 



Fig. 11. Future evolution of e -"^ ~ ui (left) and ip/{eH) ~ (T2 (right) for a modulus 
with p = 10"^, (^0 = 0.25, and U = e"^'^. 

02 ~ i72 _ ^-2^ Defining A(°) = \H-x^^'^\/H and A^ = \H-x\/H, where a; 
is given by Eq. (42), one has that A*^°) < 4 x 10~^ and A"^-*^) < 4 x 10"'' around 
z < 1. The approximated solutions deviate more than 1% from the numerical 
one when A^o) > 10"^ for z > 1.6 and A^^) > lO'^ for z > 2.5. 

As anticipated, this model-selection procedure is not sufficient by itself to 
obtain solutions which are always stable. It turns out that Stt becomes larger 
than 1 when z > 0.2 and is sub-luminal afterwards. That happens well before 
(going to the past) the point from which the numerical solution H and the 
approximated GR one x^^\ Eq. (43), begin to diverge. That at early times one 
does not get an Einsteinian evolution is confirmed by a high-redshift inspection 
of (/9, which is not constant. 

The rising of superluminal tensor modes is associated with an increase of the 
field acceleration, since Stt ~ ~'^/{H^)- We have checked numerically that, 
when Stt ^ 1, in general there is a corresponding bump in (p. Hence, one way to 
avoid superluminality is achieved by the field slowly rolling down its effective 
potential. In fact, —(p/{H(p) is nothing but the definition of the second slow- 
roll parameter ?7sr which, in this form, is independent from the type of scalar 
equation of motion. In the theory under study, (p is non-minimally coupled and, 
in general, slow rolling (r/sR < 1) will not be associated with a flat potential. 
Another possibility might be to have an effective potential with a false vacuum 
in which ip sits for sufficiently long time (i.e. with tunneling decay rate ~ H~^) 
to allow for a standard GR evolution. At some point the field tunnels down 
to the true vacuum through a potential barrier Af/ = 0{(3), and the GB term 
becomes dynamical (and, in the best case, also dynamically important). 

Both these scenarios require a complicated effective potential and a highly 
nontrivial dynamical behaviour. Without a more solid theoretical implemen- 
tation at hand, we can regard them only as proposals for future directions. One 
possibility would be to consider flux compactiflcation models [53], where the 
dynamics are given in terms of polynomial interactions between the different 
moduli (plus, possibly, nonperturbative effects). Those would give rise to more 
involved dynamics than the models presented here, and false vacua would al- 
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0.18 0.2 0.22 0.24 0.26 0.28 0.3 

Fig. 12. Region (black) in the c^o-lii/3 plane satisfying the no- ghost and sub- luminal 
conditions for the positive branch ipQ of the NOS model for f2 = U = e~^^'^. The 
contour does not change if acceleration bounds are taken into account. 

most certainly exist. It must be noted that, in its vast majority, these models 
have been constructed in the context of type IIB string theory, for which a GB 
term is not realised. However, analogous constructions have been proposed in 
the context of heterotic string [54,55,56], and it would be interesting to study 
the implication of the presence of flux in the dynamics of extended gravity 
models. 



6.2 The pseudo-modulus of the NOS model 

In the model actually considered in [13], /2 is an exponential instead of a 
hyperbolic cosine and u = 1. We have checked that even in that case there is 
a nontrivial parameter space, which is shown in Fig. 12 

6.3 Dilaton case 



On can perform the same analysis for the dilaton. We have tested a scalar 
field with 7j2 = — v^ for (i) f/ = 0, (ii) U = fi, and (iii) U = to. Since none of 
this cases is robust against theoretical and experimental bounds, as stressed 
throughout the paper and here reconfirmed, we only quote the main results 
without producing any figure. 

In general, there is a nontrivial parameter space compatible with the four 
no-ghost constraints. When imposing also the sub-luminal constraints, only 
case (iii) for the 'positive' branch (which does not give a cyclic evolution) 
survives the selection, due to the fact that Stt > 1 typically in these models. 
When increasing the magnitude of the |7i2| coefficients there appear regions 
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in which Stt < 1, but then Ssc is too large, implying that Eq. (82) is far 
from guaranteeing sub-luminal propagation. Even examples within the allowed 
regions, however, are superluminal in the past. Cases (i) and (ii), as well as 
others we have sampled, never accelerate, as anticipated in Section 3.1. 



6.4 An example of viable model 



One might ask if there is a chance to find a viable model which fits data and 
satisfies the no-ghost and sub-luminal constraints. Of course, to set /3 = 
solves the problem, but this is tantamount to saying that GB actions are not 
viable in cosmology. However, it would be interesting to find, as in the above 
analysis, an upper bound for /3 which is physically consistent. We have studied 
the evolution of these models from a chosen redshift Zj > 0. In order to do this 
we have used the number of e- foldings N = ln{l+z) as the new time variable, 
and u = ln{H/H()) and (/? as dynamical variables. As initial conditions, at 
Ni = ln(l + Zi), we have chosen Ui = \n{HciR/ Hq) and ip such that a ~ czgr, 
where the subscript 'GR' denotes the general relativistic value of the above 
quantities. d(p/dN is found by using the Friedmann constraint. In this way, the 
dynamics were initially close to pure GR (universe filled with radiation, dust, 
and a cosmological constant). For the kinematic evolution of the universe to 
be acceptable, ao > and u needs to vanish at A^ = 0. 



For the dilaton (721 = 1/8, 741 = 0.7, 742 = ±1) we found that it is, in general, 
hard to find good agreement with GR evolution, leading this to the cyclic 
behaviour already discussed before. On the other hand things change when 
we consider a modulus-like field. In fact, we have found a model which works, 
by choosing /i = 1, /2 = exp(-y3v9), 741 = 0.7,742 = 1, /3 ~ 10"^^ and (fi = 
0.125 as initial condition at Ni = 10 (zj ~ 2 x 10*^). Then the universe follows 
a GR evolution, accelerates today and all the no-ghost constraints hold. Both 
the speeds of propagation differ from one at most by one part in 10^ at high 
redshifts (see Fig. 13 for all these features). It should be noted that it is crucial 
to choose the yj^-branch as initial condition for 0, otherwise superluminal 
modes appear at late times. This means that the result is strongly dependent 
on the initial conditions. Furthermore, as (3 increases, in order to have early- 
time agreement with GR, also the value for ipi needs to increase: for example, 
for (3 = 10~^°, we require ipi ~ 4.5 at Zi, but today H diverges from Hq. 

In general, there is an upper value for /^eg, larger than the real coupling P, 
below which one recovers, for some models, the Einstein-Hilbert evolution 
from some redshift on. For other choices of the parameters we have also checked 
that, when P^s is lowered, the onset of superluminality is shifted backward in 
time and one gets closer and closer to a GR behaviour. However this means 
that, in order to get a GR evolution at early times (when the contribution of 
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Fig. 13. Comparison of the time-evolution for GR and GB cosmologies (top), a 
(central), and the speed of propagation for the scalar (dashed) and tensor (solid) 
modes (bottom) for the viable example described in the text. 

the quadratic curvature term is comparable with the linear EH), the effective 
coupling should be PcS ~ P- We finally noticed that, if /2 is assumed to be the 
hyperbolic cosine of the genuine modulus case, superluminal modes appear 
again for the same initial conditions. 



7 Including a second scalar field 



One may suggest that the approach followed so far may not be sufficient to 
select models with a runaway scalar coupled to the GB combination because. 
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in general, there will be more than one field having similar couplings. An 
example can be a complex modulus field instead of a real one. Just to give 
an idea of what changes when considering two real scalars coupled with a 
GB term, in this section we shall provide only a partial analysis of such case, 
concentrating on the behaviour of the tensor modes and giving a numerical 
example to evaluate the impact of a second field. The answer is that, at least 
when the extra field has a vanishing potential, its effect is negligible. 

We call S^ the action given by Eq. (5) with /i = 1, uj{ip) = const, and 04 = 0. 
As discussed in [1,57], one introduces the following extra terms: 

S = S^ + S^, (98) 

5^ = |rf^xv^[-i^V„xV"x + /3(x)>CGB] , (99) 

where x is another scalar modulus and ^ is constant. According to Eq. (99), 
the two scalar fields do not interact at tree level. Then the Einstein equations 
have now some additional contributions, 

+im„(^v"v,)h , (100) 

and the equation for x reads 

^nX + /3,x'CGB = 0. (101) 

In a FRW background these equations read, after redefining x ~^ v3x//t , 

f3 ih + h) H^ + H' = pm + Pr + luj^' + l ^X' + U, (102) 

ui^ + SHip) + U^ - Ph,^H\H^ + H)=0, (103) 

^(X + SHx) - (3f3,xH\H^ + H) = 0. (104) 

Since the tensor modes do not mix with scalar ones, the conditions TTl, TT2, 
and SLl can then be applied immediately, by replacing /2 with /2 + /s: 



0<l + /3(/2 + /3), (105) 

0<^±^i4±^<l. (106) 

The introduction of a new field enlarges the available parameter space, and 
a full numerical analysis is beyond the scope of this paper. However, a few 
insights in this model can be obtained by looking at the particular case 
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f2{ip) = 721 exp(722 V?) , (107) 

Mx) = 751 cosh(752 x) , (108) 

while U is given by Eq. (27). Furthermore, there are three free independent 
initial conditions, ipo, Xo, and xo, whereas 0o can be found by using the Fried- 
mann constraint (102). By choosing values of the parameters to be of order 
one (721 = 1, 722 = 2, 741 = 0.7, 742 = -2, 751 = 752 = 1) and initial condi- 
tions <^o = Xo = 0, the field x evolves to a constant value, i.e. x ~^ 0. This 
solution is an attractor, because it is reached independently from the initial 
condition given to Xo- Therefore, at late times the evolution tends to be the 
same solution of the NOS case discussed before. In particular, as it happened 
in the NOS case, the tensor modes are superluminal. 



8 Discussion and conclusions 



From the analysis of different existing models we conclude that the appearance 
of ghosts and other particle instabilities, as well as of causality and Lorentz- 
invariance violation, is not an uncommon feature in Gauss-Bonnet gravity. 
In particular, we have seen that few of the string-inspired models examined 
so far are theoretically viable. It is still unclear if there is a nonempty set 
of initial conditions satisfying all the no-ghost and sub-luminal conditions 
during the entire evolution of the universe, but it has been established that 
these constraints do limit, in a decisive way, the parameter space of initial 
conditions. 

Cosmological solutions which are phase-space attractors are not guaranteed to 
be ghost free. Stability in the phase space {H, H, 0, 0} against homogeneous 
perturbations does not imply unitarity at the quantum level. Our analysis 
covers cosmologically long time intervals, and the possibility of allowing for 
the existence of unobservable, decaying ghost modes is not realistic within our 
framework. This fact contributes to ruling out most of the models studied. 
Some of these are indeed ghost free at any time but perturbations propagate 
at a speed faster than light. We have also given an example of a model viable 
from z ~ 10*^, but it is motivated only phenomenologically. 

Therefore, we can conclude that the presence of ghost modes for the dilaton 
and moduli fields may be a compelling reason to stabilise them before they 
can play some role in cosmology. Our results indicate that the non-minimally 
coupled scalar field of the action (5) should acquire a precise value through 
the cosmological dynamics. The alternative to face is to carefully fine tune the 
initial conditions of the equations of motion, which makes this class of models 
rather unattractive. However, its study is far from being complete. 
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The presence of ghosts and superluminal modes in effective Lagrangians ob- 
tained by string theory on FRW backgrounds should not surprise the reader. 
In fact, all these models are based on an effective Lagrangian which was found 
only assuming a Minkowski background. However, as we have shown in this pa- 
per, the same Lagrangians cannot be used tout court in cosmology, otherwise 
unphysical modes may propagate. On the other hand, an effective Lagrangian 
should not depend in principle on the particular background we choose to 
study: even the Einstein-Hilbert action is used to describe any kind of gravi- 
tational system. It is then possible that cosmology may actually help exclude 
those actions which give rise to unphysical states. 

To achieve quantum stability and a satisfactory cosmological evolution it is 
necessary to impose a particular form for the potential. In our case we must ei- 
ther require some nonperturbative mechanism that produces potentials which 
can heal the FRW (and any other) background from ghosts and superluminal 
propagation, or we should regard the GB action as an effective one to be cor- 
rected by loop contributions. These were found for type II string with fixed 
moduli in [58,59,60,61] and in dilaton super and heterotic string in [21,62]. The 
results are summarized in [63]. We note however that, in general, there are 
terms which add to higher-order Euler densities and lead to effective theories 
with spin-2 ghosts even on Minkowski background. As stressed in the introduc- 
tion, higher-order curvature terms are subdominant at low energies anyway, 
while at high energies a new and more involved analysis, both of background 
and perturbed equations from a nonperturbative action, is required. 

We have not fully addressed the issue of the acceleration of the universe. 
Equation (54) gives the necessary information for constraining the evolution of 
a and so that to reproduce the present inflationary expansion. Although we 
have performed an analysis of such dynamics in comparison with cosmological 
data, it is worth further investigation. In particular, we were able to flnd a 
ghost-free and sub-luminal and today-accelerating solution, but not to fully 
constrain the parameter space to be compatible with experimental bounds 
from nucleosynthesis until today. We have considerably restricted it, though. 

In order to study more realistic (from the string theory point of view) scenar- 
ios, one should extend the couplings appearing in Eq. (5) to be dependent on 
more than one (complex) moduli. This has recently been considered in [57], 
and here we have addressed it in a very simplifled way. The main conclusion 
is that the ghost constraints from the tensor modes are unaltered, and those 
are enough to rule most models out. However, it would be interesting to study 
more complex scenarios with several moduli fields evolving simultaneously. 
This immediately leads us to an interesting future direction to explore which 
is to study how scalar potentials more complicated than Eq. (27) affect the 
dynamics of the system. From our analysis there is evidence that exponential 
potentials are hardly able to describe cosmologies which are both free from 
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instabilities and compatible with observations and, in general, one is forced to 
trade one shortcoming for the other. Considering flux compactification mod- 
els would combine both the inclusion of more than one scalar field into the 
analysis, with more involved dynamics parametrised by the flux terms in the 
scalar potential, and would certainly lead to interesting situations. 

We have chosen a specific form for the functions /», but these can be more 
general beyond tree level. For example, one may invoke exponential series from 
higher- loop corrections, as in [30], or nonperturbative effects. It is clear that 
the viability of these models depend on the choice for /j, and a number of 
cases of interest is still to be explored. 

Sub-luminal propagation of perturbative modes is one of the conditions which 
limit the parameter space of the theory. In scenarios where faster-than-light 
travel is allowed, such as theories defined on non-commutative backgrounds 
[64] , we would expect the space of initial conditions to be modified and some- 
how enlarged. However, the parameter space would be still non-trivial. This 
is also true when, for some other reason, one abandons the sub-luminal con- 
straints altogether. 

Finally, if matter and radiation are non-minimally coupled with the cosmo- 
logical evolution may change considerably [14]. 

In this paper we have established, both analytically and numerically, the fol- 
lowing very general results: 

• Classical and quantum stability are inequivalent issues which require sepa- 
rate assessments, one not implying the other. Most of the models presented 
in literature were tested in phase space, but for arbitrary initial conditions 
they have ghosts or other kinds of particle instabilities, as well as problems 
related to the propagation speed of particle degrees of freedom. 

• In particular, even the most orthodox models based on string theory (i.e., 
flat directions like the dilaton or a compactification modulus) are not auto- 
matically ghost-free: we are forced to select special initial conditions. 

• It is very difficult, although not impossible, to obtain a theoretically and 
experimentally viable picture from this class of models. We were able to give 
only one positive example, with a GB coupling apparently well below any 
reasonable threshold of experimental detectability. Still, minimally coupled 
General Relativity is by far the most successful theory at our disposal. 

• Qualitatively, the presence of many noninteracting scalar degrees of freedom 
does not change the main features of the theory, including the existence of 
nontrivial selection rules on the parameter space. 

The presented parameter selection procedure, although remarkably effective, 
clearly falls short of being fully prescriptive, and there is some evidence that 
non-minimally coupled GB cosmology is unlikely to be a viable theory for 
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natural values of the parameters and simple choices for the scalar functions. 
We hope to further clarify this and at least part of the above points in the 
near future. 
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